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(-H . Abstract 

. _ The propagator for a certain class of two time-dependent coupled and 

■ driven harmònic oscillators with time-varying angular freqüències and masses 
^ . is evaluated by path integration. This is simply done through suitably chosen 

generalized canonical transformations and without presupposing the known- 
ledge of any auxiliary equation. The time-dependent oscillators system with 
^ ' an exponentially growing masses and coupling coefficient in time may be 

■ considered as particular case. 
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PACS 03.65. Db - Functional analytical method 



1 



^The Feynman path integral formalism is known to be an elegant and 
powerful approach for analytically treating some time-dependent one dimen- 
sional physical systems. We can quote , for instance, the harmònic oscillator 
with variable frequency or mass and a perturbative forcdïliffli , the infinite 
potential well with moving wall,i the time-dependent potentialsB of particular 
type V{x — f{t)) and the time-dependent modeli of the form V(g/ç(í))/ç^(t) 
with ç(í) = (aí^ + bt + cY^"^. We wish to add to this list the system of two 
time-dependent coupled and driven harmònic oscillators. No exact solution 
of this problem by path integration is known to us. On the other hand, there 
has been two attempts to solve this kind of problems. The first one has been 
devoted to a system composed by two time-dependent coupled oscillators 
with damping and time-dependent external forcei . The resolution of the 
Schròdinger equation has been made by means of the construction of a time- 
dependent invariant which is linear in position and momentum variables, but 
the solution has not proved satisfactory since, the coupling coefíicient does 
not appear in the wave function. The other one relative to a system of two 
time-dependent coupled and driven harmònic oscillators with time-varying 
angular freqüències and masse . The propagator has been built through a 
variant of the sm(1, 1) Lie algebraic approach, but it is only vàlid in the case 
where the masses are constant and starting with a rotational transformation 
parametrized by a constant angle. 

Our intention in this paper is to calculate the propagator for this system 
within the framework of the path integral approach. The calculation is facil- 
itated by utilizing the Hamiltonian formalism and the generalized canonical 
transformat ions . 

For the dynamical system of our interest (see Ref.0), the Hamiltonian 
has the form: 



p^- 1 



+ \{t)xiX2, (1) 



where mj(t), ujj(t), fj{t) and X(t) are the time-dependent mass, angular fre- 
quency, force and coupling, respectively. This Hamiltonian is a generalization 
of the Hamiltonian system considered in Ref.i , where trijit) = rrije'* and 
\{t) = rmim2e^^. 

The quantum mechanical evolution of the system can be described by 
the propagator, in the phase space formulation of Feynman's path integral. 
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which is defined formally by 

K{x'l, X2, t"; x[, X2, f) = J VxiVx2VpiVp2 exp j^sj , (2) 

where 

S = J^, ÍPi X, +P2 X2 -H{t)) dt (3) 

is the Hamilton's principal action integral. 

In a time-graded representation, this expression (0) is understood as: 



Af-l N 



K{x'l,X2,t";x[,x'2,t') = lim / TT dxindx2n n 



dpin dp; 



'2n 



n=l 



^^^P\j-Y.Sn\, (4) 



n=l 



with a short-time action 

•S'n = (pinAxin + P2n^X2n " sHiJn)) ■ (5) 

Here, we have adopted the Standard notation: e = t„ — í„_i = {t" — t') /N, 

tn (tn "t~ ^n— l) /2, 2;^^ Xj(tn): Xj Xj(t ) Xjitç)), Xj Xj(t ) 

a^iíií^Tv); i = 1,2. 

The path integral is not trivial but can be evaluated. In order to cast 
the integral into a more tractable form, it is advantageous to use the method 
of canonical transformations by performing the time-dependent canonical 
transformation {xi,X2,Pi,P2,t) — > {Qi,Q2, Pi, P2,t) defined as: 

Qi cos a{t)+Q2 sin a(t) — Qi sina(t)+Q2 cosa(t) 

pi = Jmi{t) (Picosa(í) + P2sina(t) + p^{t)xi) , (6) 



P2 = ym2{t) {-Pi sina(í) + P2COSQ;(t) + (32{t)x2) , 

where the functions and [32{t) can be conveniently chosen to make 

separation of variables straightforward possible. 
From the classical mechanics equationJHI 

Pj = ^F2(xi,X2,Pi,P2,t), = ^P2(a;i,X2,Pi,P2,t), 3 = 1,2, (7) 



and 



n(t) = H(t) + |f2, 



(8) 



the generating function responsible for the transformation is found to be 
F2{xi,X2,Pi,P2,t) = Jmi{t) {Pi cos a{t) + P2sma{t)) xi 



+\/m2{t) (— -Pi sin a{t) + P2 cos a{t)) X2 
I3^{t)^jmi{t)xl + P2Ít)\lm2{t)x 



1 

+2 



(9) 



In the new conjugate variables {Pi, P2,Qi,Q2), the Hamiltonian for the 
system becomes 



n{t) = i (P^ + P|) + A{t)P,Qi + B{t)P2Q2 + {P,Q2 + P2Q1) 
+ « (í) - + \Di{t)Ql + ^L'aWQa + E{t)QiQ2 



-Fi{t)Qi - F2{t)Q2, 



(10) 



where the time-dependent coeficients A{t), B{t), C{t), Di{t), D2{t), E{t), 
Fi{t) and F2{t) are given as 



A{t) 



COS^ «(í) + 



^mi(t) 2mi(í)y \^^m2(í) 2m2(í)^ 

= + "^iW ^ sin^ aíí) + í ''^aft) , "Míl^ cos^ a(t) 



1 / rni(t) , m2(t) \ , / /3i(t) _ feW 

2 V"ii(í) m2(t)y \^mi{t) s/miit) 



sinQ!(í) cosa;(í) 



Di{t) = di(í) cos^ «(í) + d2{t) sin^ «(í) - 



m2(í)^ 

2 „yj.\ 2A(t) sina(t) cosa(f) 



D2{t) = di{t) sin^ «(í) + d2{t) cos^ «(í) + 



mi{t)rn2(t) 
2A(í) sin a(í) cos q(í) 



E{t) = {di{t) -d2{t))sma{t) cosa(í) + 



A(í)(cos2 Q;(t)-sin2 a(í)) 
^/mi{t)m2 (í) 



Fi(í) = ^Jm^fl{t) COS a{t) - ^Jm^f2{t) sin a{t) , 
F2{t) = ^mi{t)h{t) sina{t) + y^^/2(í) cosq;(í), 



with 



mi(í) mi(í) 
4 



mi(í)/3i(í) 



(11) 



(12) 



and 



d2{t) = ul{t) + 



+ 



1 d 



m2(t) Tn2{t) dt 
As a result of these transformations, the action (Q) takes the form: 



(13) 



S 



t" 



-Qi Pi -Q2 P2 -n{t) + 



dF2 
dt 



dt 



(3^{t)^mi{t)xl + (32{t)^m2{t)x 
+ £ [p,Q,+P2Q2-n{t))dt. 



(14) 



By means of a procedure similar to that presented in Ref.i , the measure of 
the path integral in (H) changes as 



r 1 - 

VxiVx2VpiVp2 =n Hi kJ ' T^QpPj^ 



(15) 



where m'- = mj{t") and m'^ = mj{t'). It therefore follows from (0) and (15) 
that the propagator (^ may be written in the form: 



K{x'l,x'i,t";x[,x'2,t') = n 



exp 



2h 



^^{t)^m,{t)x] 



xK(Q^,Q2,t ;Q^,Q2,t^ 



(16) 



where 



K{Q'U Q'^, t"; Q[, Q'2, t') = Jll VQpP^ exp |^ (Pi +P2 Q2 -n{t)) díj 

(17) 



is the propagator for the physical system governed by the Hamiltonian 7i(t) 
in the new conjugate variables. 

If we make the foUowing choices: 



a(t) = Const, 



(18) 
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and 



/3,(t) 



(19) 



the terms in P1Q2 and P2Q1 are cancelled and the new Hamiltonian ([ÏOD 
takes the form: 



n{t) =E 



+ r(í)QiQ2, 



(20) 



where 



íli{t) = iül{t) cos^ a + iül{t) sin^ a 7 

Q2Ít) = J^\{t) úv? a + ülit) cos^ a + 



A(í) 



mi(t)m2(í) 



sin(2Q;), 



A{t) 



-y^mi (í)m2 (í) 



sin(2Qí) 



Fi(í) = mi(t) fi(t) cosa - ^Jm^f2{t) sina, 
F2(t) = ^mi{t)fi{t) sina + \Jm^f2{t) cosa, 



(21) 



r(í) = i (í) - íl)^(í) sin(2a) + 



A{í) 



Y^mi(í)m2{í) 



with 



C?(í) 



cos(2a), 



J = l,2. 



(22) 



4 l rrí^it) rrijit) 

Notice that, with the above canonical transformation, the couphng T{t) is a 
functional on the parameters of the original system. Then, the separation of 
variables occurs for a wide but restricted class of two coupled general driven 
time-dependent harmònic oscillators contrary to this what has been asserted 
in Ref B . 

It is clear that the separation of variables in Eq. (pÜ[) requires that V{t) = 
0, i.e. 

A(í) = ^^mi{t)m2{t) (ul{t) - ül{t)) tan(2a). (23) 
Then, by taking into account (p3D, the propagator (|T^ is rewritten as: 

K{Q';,Q'^,t";Q'i,Q2,t') =n K{Q';,t";Q'j,t'), (24) 



where 



K{Q';,t";Q'^,t') 



VQpP, exp{- Pj Q. -Hj (t) dt 



(25) 



is the propagator for a time-dependent harmònic oscillator with a perturba- 
tive force described by the Hamiltonian 



Pf 



H,{t) = ^ + -ü'^{t)Qj~F,{t)Q,. 



(26) 



Now, it is easy to write down the propagator (^). Using the generahzed 
canonical transformations method, Chetouani and co-workerJÍ have devel- 
opped a very elegant way of calculating the propagators of certain time- 
dependent one-dimensional physical systems. To use this method, we apply 
the canonical tansformation, followed by a time transformation, defined by 



ds 
dt 



Pfit), 



(27) 



where Pj(í) is an arbitrary function without dimensions. After some calcula- 
tion, we find that 



K{Q''t";Q',t') 



where 



1 




K{Xls"-X'^,s') 



VVpV, exp - 



p2 

3 



^3 - hf 



%{s)p\\X'^-F,{s)-p]XAds 



(29) 



with 



Pi 



^ -2 ( ^ 

P, 



P] Pj 



In the above we have used the notation = Pj{s) 



FJs) 



Fj{t) and í]2(t) =.n'^.(s). 



Pj= Pf 



(30) 



ds ' 



A comparaison between the propagator (^) and (^) shows that the 
space-time transformations (0) has resulted in the appearance of a phase 
and a quadratic term ^Ú'jXj. On the other hand the global time-dependent 



frequency in Eq. 



depends on the adjustable parameter p . 
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If we now impose a constraint on by setting the global time- dependent 
frequency appearing in Eq. (pÜ]) equal to a constant: 

+ ü] {s)fj = u% = Const, (31) 



we recognize the propagator of Eq. (|29|) as the expression given by Feynman 
and HibbJi for the forced harmònic oscillator with a constant frequency and 



we therefore find for the propagator (p8 



• // 



X exp 



2hsm(l)j{t",t'' 
PjPj ^Oj Pj Jt' 

+ — / GAt)s\Ti(t)At",t)dt 

UJnn Da Jt' ■' 



UjQj Pj 

2 /•*" /•* 



ÜÜQj Jt' Jt 



Gj{t)Gj{T) sin ^At",t) sm(j)AT,t')dTdt 



(32) 



where 



G,{t) = F,{t)p^{t) and 0^.(t",í') = o^o, j^^ - (33) 



it' 

Substituting ( PUD into (0), we see that Pj(t) is solution to 



h +^mp. = (34) 

Pj 

which is the well-known auxiliary equationSB . 

Inserting (^) and (^4]) in (fLÜj ) and using (p^Üf), the final expression for the 
system of two time-coupled and driven harmònic oscillators governed by the 
Hamiltonian is given by 



iÜQj 



\ \ 2innp]p'jsm(f)j{t",t'] 



exp < — — 



ih rrijit) ' 
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/ ■ II I 

Pi ^119. P 



X 



exp 



j n"2 o /n'2 l 



2^sm0.(t",í') 



— ^ + -^ COS0^(í",í') 

Pj Pj / 



2g''g' 2 Q" rt" 



+ 



Gj(í) sin 



+ 



2 Q' 



UjQj Pj Jt' 

2 /-í" rt 



Gjit) sin (f)j{t",t)dt 



— 2" / / Gj(t)Gj{T) sin (f)j{t" ,t) sin (j)j{T,t') drdt 

UÍQj Jt' Jt' 



(35) 

Next, we turn our attention to the special case where mj{t) = rrije^^ and 
A(t) = rmim2e"'*. The propagator defined by expression (0) can be written 

K{xl t"; x[, x'„ t') =n /C(Q;', t"- Q;, í'), (36) 

i=i 

with 



K:(g;',í";g;,t') 



2mha'-a'jSÍn(t)j{t",t 



X 



exp 



2;isin0j.(í",í') 



Q 



,112 



+ 



cos 0,·(t", í'^ 



+ 



2 



2 q;- 



Gj{t) sin ^At,t')dt 



Gj{t) sin (j)j{t",t)dt 



2 /■*" /■* 

— 2- / / Gj(t)Gj{T) sin(f)j{t" ,t) sin(j)j{T,t')dTdt 

LÜQj Jt' Jt' 



(37) 



where aj = pj / ^JmTy 

In the absence of the couphng, i. e. A(í) = and furthermore a = 0. 
Our system is then composed by two uncoupled harmònic oscillators with 
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time-dependent and different angular freqüències, masses and forces. In this 
case, the propagator is 



K[xi^, X2,t ; X]^, X2, t ) —Y\_ -^i-^j ' ^ i -^j' ^ ) 



(38) 



where 

K{x';,t";x'^,t') 



iüoj 



2inha'^a'jSm(l)j{t",t' 

iuJoj 



cr 



3 



Oi 



X exp 



^ + -12] COS <f)Jt",t') 



a 



a 



2t"-t'- 9 r" rt" 

+ I G,(t)sin0,(M')cit 

-—^t GAt)ún<pAt\t)dt 

üjQj (Tj Jt' 

- — / / G'j(í)G'j(r)sin</)j·(t",í)sin</)j·(r,í')c?^c?í 

Cc'gj Jt' Jt' 



(39) 



is in agreement with the one obtained in Ref.i . This is clearly proof that 
our exact results ( P^ and (^) are indeed the correct quantum propagators 
for two time-dependent coupled and driven harmònic oscillators. 

In conclusion, this presentation contains two principal results. First, we 
have succeeded in extending the path integral method to include a wide class 
of a pair of time-dependent coupled and driven harmònic oscillators. By 
explícit path integration, we were able to derive closed form expressions for 
the Feynman propagator in a rigorous and explícit way. 

Second, the correct propagators for the systems included in this study 
have been obtained for the first time. 
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